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Abstract 

The extension of the Lieb-Schultz-Mattis theorem to dimensions larger than one 
is discussed. It is explained why the variational wave- function built by the previous 
authors is of no help to prove the theorem in dimension larger than one. The short 
range R.V.B. picture of Sutherland, Rokhsar and Kivelson, Read and Chakraborty 
gives a strong support to the assertion that the theorem is indeed valid in any dimen- 
sion. Some illustrations of the general ideas are displayed on exact spectra. PACS 
numbers: 71.10.Fd; 75.10. Jm; 75.40.-s; 75.50.Ee; 75.60.Ej; 75.70.Ak 

1 Introduction 

In 1961 Lieb, Schultz and Mattis proved that a spin 1/2 antiferromagnetic peri- 
odic chain of length L has a low energy excitation of order 0{\/L) Q. This theo- 
rem (called in the following LSMA) was then extended by Lieb and Affleck to odd- 
integer spin but fails for integer ones Pj. It states that SU(2) invariant Hamiltonians 
with odd-integer spins in the unit cell, either have gapless excitations or degenerate 
ground-states in the thermodynamic limit. The authors suggested that it might be 
extended to higher space dimensions, but up to now, no complete argument has been 
worked out ||. 

In this paper we revisit the method used by LSMA (construction of a variational 
excited state) and the physical meaning of the unitary operator involved in this 
construction. This sheds some light on the reason why the LSMA excited state is 
generally not a low energy excitation in space dimension larger than one (as for 
example in the case of Neel order on the triangular lattice), and how qualitatively 
one might try to build a truly low energy excited state. We then study an alternative 
wave-function very much in the spirit of the resonating valence-bond (RVB) states 
of Sutherland ||, Rokhsar and Kivelson ||, and Read and Chakraborty ||. In this 
last framework, we show that the above-mentioned statement is indeed true in any 
dimension, and discuss the quantum numbers of these quasi degenerate ground-states. 
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2 The LSMA Theorem 



To analyze LSMA theorem in dimension larger than one a dichotomy could be done 
between situations with T = long-range order and symmetry breaking order param- 
eters on one hand, and systems without long-range order on the other. 

• In the first case the theorem is trivially true: a symmetry-breaking situation 
necessarily implies in the thermodynamic limit the mixing of states belonging to 
different irreducible representations of the Hamiltonian and / or lattice symmetry 
group and thus a degeneracy of the ground-state. 

i) In the case of Neel long-range order, there are both degenerate states (which 
form the true thermodynamic ground-state) and gapless excitations. The gap- 
less excitations, the antiferromagnetic magnons, are the Goldstone mode asso- 
ciated to the broken continuous SU(2) symmetry. Whereas it is well known that 
the softest magnons scale as 0(1/L), it is often asserted that these quasi-particle 
first excitations are the first excited levels of the multi-particle spectra. This is 
indeed false: the "T=0 Neel ground-state" (or "vacuum of excitations" ) is itself 
a linear superposition of ~ N a eigen-states of the SU{2) invariant hamiltonian 
which collapse to the ground-state as 0(1/ L 2 ) (a is the number of sub-lattices 
of the Neel state, L is the linear size of the sample and N the number of sites 
of the sample) 0, |, @. 

ii ) In the case of a discrete broken symmetry, as for example the space symmetry 
breaking associated to long-range dimer or plaquette order, or in the case of 
T-symmetry breaking associated to long-range order in chirality, the vacuum of 
excitations is also degenerate in the thermodynamic limit. The number of quasi 
degenerate states is related to the dimension of the broken symmetry group, 
the nature of the order parameter and is independent of the sample size [|lC|], 
The collapse of these levels on the ground-state is supposed to be exponential 
with the size of the lattice. 

• It is only when there is no long-range order, a spin gap and when all spin 
correlations are short-ranged that the existence of a ground-state degeneracy 
is non-trivial. In that case, the theorem states that the state(s) degenerate 
with the absolute ground-state in the thermodynamic limit have wave vector 
7r with respect to the ground-state. More generally, in dimension larger than 
1, we will argue that the states which collapse to the ground-state have wave 
vector (wave vectors joining the center to the middle of the sides of the 
Brillouin zone) Q. These low lying levels collapse to the absolute ground-state 
exponentially with L. 

3 The LSMA variational state 



1 This result supposes that the ground-state wave vector is zero. It is zero in all situations we have been 
looking at. If a system had a non-zero momentum ground-state, it would be degenerate and the theorem 
true anyway ! 
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3.1 Finite size spectra and Twisted Boundary Condi- 
tions. 

The spectrum of a given Hamiltonian on a finite sample depends on the boundary 
conditions: these can be free, periodic or twisted. The LSMA theorem is valid for 
periodic boundary conditions. For a complete understanding of the nature of the 
LSMA variational state it is useful to study how a given exact spectrum evolves 
under a twist of the boundary conditions. 

Let us consider a finite sample in d dimensions, described by d vectors Tj. Gen- 
eralized twisted boundary conditions are defined by the choice of the twist axis (here 
the z axis in the original Bq spin frame) and a set of d angles <f>j as: 

S(Ri + Tj) = e ^ sz ( R »)s(Ri)e-^ Sz(Rl) (1) 

From now on, and for the sake of simplicity, we will develop the algebra on simple 
1-dimensional models (extrapolation to larger dimensions or more complicated SU(2) 
invariant Hamiltonians is just a problem of notations). As an example, let us consider 
the nearest neighbor Heisenberg Hamiltonian with periodic boundary conditions: 

L-l 

7~to = ^2 Sn ■ S n+ l[L] (2) 

71=0 

A twist 4> hi the boundary conditions implies the calculation of the eigenstates of 

L-2 

= Si_i ■ (e^Soe-^o) + £ S n • S n+1 

n=0 

= n + l((e»- 1)^-^0+ +h.c.) (3) 

Under an adiabatic twist of the boundary conditions the spectrum of Ti^ evolves 
periodically with a period 2tt, as the boundary conditions (Eq. |l|). But the eigenstates 
evolution might be more complicated: a unique spin-^ wave function acquires a phase 
factor —1 under a 2tt twist. And there is no guaranty that the the ground-state of 
7~C</)=o adiabatically transforms into the ground-state of H$=2-k- As we will show 
below, this is generally not the case and the true period of the eigenstates is 4-7T. 

3.2 Twisted Boundary Conditions and Translational In- 
variance 

To follow adiabatically an eigenstate while twisting the boundary conditions may 
be difficult if there are level crossings during the twist. The only way to do it in an 
unambiguous way is to follow a given one-dimensional representation of the symmetry 
group during the twist: in such a representation the levels are non degenerate and 
never cross. The ground-state of the Hamiltonian with periodic boundary conditions 
indeed belongs to such a representation 0. 

2 Up to now, in any antifcrromagnetic models we have studied, the ground-sate of an even number 
of spins belongs to the trivial representation of the total group (the only exceptions are associated to 
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Such a program is not easy in the framework of Eq. ^, which breaks the transla- 
tional symmetry of the problem. But the translational invariance can be restored by 
a unitary transformation rotating the spin frame at each lattice site. Let us call Ba,, 
the new frame deduced from the original Bo by a spatially dependent twist described 
by the unitary operator: 



U(4>) = exp(i-^ nSf, 



n=0 



In this new frame, the twisted Heisenberg Hamiltonian reads: 



(4) 



(5) 



(in this equation and in the following, we put a tilde on each quantity measured in 
the Ba, frame). This unitary transformation is chosen so that the boundary term(s) 
in Eq. || disappear (s): 



L-l r 

£ 

n=0 



S n S n+l[L] + 2 [ e% ' t>/LS n S n+1[L] + h - C - 



7~Cd> = 7~Lq + 



,L-1 
n=0 



[(e^/ £ 



+1[L] 



+ h.c. 



(6) 



(7) 



Ti.^ is translation invariant (|?1^,T = 0, where T is the operator for one-step trans- 
lations) and its spectrum is indeed identical to the spectrum of TIa,. We can now 
define irreducible representations of the translation group labelled by their wave vec- 
tors in the Ba, frame and adiabatically follow a given eigenstate of the momentum in 
the successive Ba, frames while increasing <p from to 4-7T (see example in Fig. ||). 



For a given twist 4>, the zero-momentum eigenstate 
has for expression in the Bo frame: 

i>f) = U-\cj>) $f 



k=0 



of Hri, in the £u frame 



(8) 



For an arbitrary twist, ip^J does not describe a spatially homogeneous state in the 
Bo frame. 

We will now show that for a 2ir twist, the trivial representation of the translations 
in the B^ frame ( ^zk/)) ^ as momentum k^. in the Bo frame. Following Affleck Q, 
this is easily done by noting that for odd-integer spins U(2ir) anti-commutes with T, 
as soon as the number of rows in the transverse direction is an odd integer [3[|. This 
proves that jV'lD defined by Eq. || takes a phase factor —1 in one-step translation 
along the twist direction and thus has a momentum k^ in the Bo frame. 



special pathological behavior of very small samples): it is a state with total spin zero, zero momentum, 
invariant in any operation of the point symmetry group. This makes sense and appears as a powerful 
extension of Marshall theorem for bipartite lattices. If the ground-state were to belong systematically to 
a multidimensional representation, as it might be the case for chiral spin liquids, then the LSM theorem 
would again be trivially true. 
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An example of such an adiabatic continuation is shown in Fig. [2[ The spectrum 
of the multi-spin exchange Hamiltonian on a small losange (4 x 5) is displayed as 
a function of the twist angle. The boundary conditions are twisted in the direction 
of length L = 4 (the number of rows is odd). We can note the above-mentioned 
properties: 

• The spectrum is periodic in <f> of period 2ir. 

• The eigenstates of and T evolve with a period Att. 

• For a 2-7T twist, the zero momentum eigenstate of H.2n hi the Bi-k frame has a 
momentum in the Bo frame (compare the spectra and labels for twists 
and 27r) 



3.3 The LSMA variational state revisited 

The proof of LSMA theorem relies on the construction of a low lying excited state for 
the problem with periodic boundary conditions. Let us call \ipo) the exact ground- 
state of this problem. The LSMA excited state is obtained by the action of the 
unitary operator U(2ir) (Eq. |j) on |V'o) : 

4 5AM ) = U(2n) |V> ) (9) 

As already mentioned, in specific geometries, U(2tt) anti-commutes with T and 
^2if Mj4 ) nas momentum it with respect to \4>o). Contrary to the states described 

above, #2t? MA ) I s n °t an eigenstate of Hr>. Its variational energy can be computed 
with elementary algebra as: 

(ei SMA \ n \ei SMA ) = m . (10) 

Using Eq. ^ (generalized to d dimensions) we see that 02^ MA ^j variational energy 
is equal to the exact ground-state energy (-00 1 ~Ho IV'o); plus a correction of the order 
of 0(N x ( e 2i7r / L - 1)), that is 0{N/L 2 ). This gives, in one dimension, an energy 
0(1/L) and achieves the proof of the LSM theorem. The number of correction terms 
does not allow to extend the proof to higher space dimensions []. 

A simple counter-example showing that d^^j is generally not a low lying exci- 
tation is given by the spectra of the Heisenberg hamiltonian on the triangular lattice: 
in this situation the exact lowest excitation energy (and thus the LSMA variational 
energy) in the sector is indeed 0(1): it is associated to the corresponding magnon 
which has a large energy in this three-sublattice ordered structure ||. 

The point of view developed in the previous subsection provides another way to 



look at the right-hand side of Eq. [L0[ [-0o) might be seen as a variational guess to solve 
the twisted boundary conditions problem using the frame, and then equivalently 
&2n MA ) as the solution of the same physical problem in the Bq frame (7^=2^ = ^o)- 
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Only systems with aspect ratios going to zero as the size goes to infinity have degenerate states 



reminding of the one-dimensional problem 
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This new point of view enlights the weakness of the LSMA wave function (in d > 1) 
and how we might try to overcome it. 

Let us for the moment assume that <j) is equal to 2ir — e. In the Bq frame, 
the perturbation to the periodic problem only appears at a boundary of dimension 
(d-1): but the variational solution ^It?-^) smears out the spin response on the 
entire system. This would be a sensible solution in the case of long range order, 
where the system shows stiffness and sensitivity to the boundary conditions. On the 
contrary, in the present case, where spin-spin correlations have a finite range £, it 
seems reasonable to expect that the boundary perturbation does not propagate at a 
distance much larger than £ from the boundary. The LSMA solution is thus probably 
very far from optimal. 

We might thus expect to find in the Bq frame, a solution with energy lower than 
qLSMA^ ky perturbing \ipo) only locally at the boundary. For such a wave function, 
Equation (|3|) (and its generalization to d dimensions) implies a distance in energy 
from the ground-state of the order of e 2 £iv We could then speculate that this 
difference in energy might be tuned to zero by an appropriate choice of the small free 
parameter e. But indeed such a reasoning involves a difficult and may be pathological 
limit. In the following paragraph we thus propose a new variational wave function 
for the low lying k,4 excited state: in the translationally broken picture {Bq), this 
excited state differs very little from j^o) an d only in the vicinity of the boundary 
defect. 

As an illustration of the present analysis, we can look at the evolution of the low 
lying levels of the spectrum of the MSE hamiltonian under a twist <j) of the boundary 
conditions (Fig. |2|). In spite of the very small size of the sample (4x5), one clearly 
sees that the exact ground-state energy does not increase as <fi 2 but more or less as 
</> 4 (see Fig. ||). For small enough <p , this state does not present stiffness to twisted 
boundary conditions: this is exactly what is expected from a Spin-Liquid. 

For comparison one can compute the variational energy of the state which inter- 
polates between |^o) and 



)LSMA\ j t ig defined ag . 



\6 4> ) = U-\<P)\^q). (11) 

Its variational energy can be rewritten as a linear combination of 2-body and 4-body 
observables of |^o) multiplied by cosines of <j>/L. It thus increases as (j) 2 and has a 
non zero stiffness, which explains why it is a bad approximation of the exact state. 

Remark: Contrary to the the TV = 36 spectrum of Fig. |l| (where £ > L), the sample 
in Fig. ^ does not display all the features expected from a Spin-Liquid spectrum. The 
correlations at a distance 4 are still not completely negligeable: this explains the rapid 
increase in the ground-state energy of TL^ for <f> > 0.2. As a consequence, the ordering 
of the eigen- levels of Hq is different from the thermodynamic limit: in particular the 
ground-state in the k^4 sector is not the first excited state of the 4x5 spectrum, 
as it is in the 6x6 example of Fig. [l|. The evolution under a twist of the N = 36 
spectrum, would have been much more pedagogic: it is for the moment too expensive 
in computer time. 

At that point it is interesting to discuss Oshikawa's approach of this question Q| . 
In place of twisting the boundary conditions, Oshikawa suppose that the quasi par- 
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tides of the problem have a fictitious charge. He inserts a magnetic flux in the hole 
of the torus defined by the boundary conditions (see Fig. [|]) and then increases the 
flux from to 2ir. This is a procedure absolutely similar to our twist of the boundary 
conditions in the spin problem, and the adiabatic insertion of a 2tt flux brings the 
system from the l^o) state to the ^=0^) °f ^2tt defined in subsection 3.2. At that 
point he does the implicit hypothesis that the gap of the system does not close in the 
operation. He thus automatically arrives to the conclusion that the ground-state is 
degenerate. 

Our approach gives a complementary insight in the physics of such a system: a 
Spin-Liquid, with its absence of stiffness is not sensitive to a twist of the boundary 
conditions (in the thermodynamic limit). Thus the ground-state energy of TL^ in the 
translation invariant sector does not depend on the twist eft, which implies that the 
ground-state of TCq is degenerate (with wave vectors and in the Bq frame). In 
the following paragraph, using a specific mathematical definition of the Spin-Liquid 
state, we will exhibit a variational wave- function giving a strong support to the LSMA 
conjecture. 

4 The short range RVB picture of the first ex- 
cited states 

In this part we use the main ideas of Sutherland [0], Read and Chakraborty 
to build an explicit variational wave-function orthogonal to the ground-state and 
collapsing to it in the thermodynamic limit. We then show that these first excited 
states have momentum with respect to the ground-state. 

Sutherland first showed that the zero-temperature observables of a nearest neigh- 
bor resonating valence bond wave-function can be computed thanks to the classical 
properties of a gas of loops. In the quantum problem, the loops appear formally 
when scalar product of wave- functions (or matrix elements of spin permutations) 
are written in terms of dimer coverings. The Sutherland nearest neighbor resonating 
valence-bond (NNRVB) wave-function description can be mapped to the high tem- 
perature disordered phase of the classical loop model: its correlation length is finite 
and the weight of long loops is exponentially decreasing with their length. 

Our own reasoning rests on the following basic assumption (A): The ground- 
state |V>o) is a R.V.B. state, and the long loops weight in the norm (^o| ipo) decreases 
exponentially with the loop length. This last requirement implies the exponential de- 
crease of all multi-point correlations with distance, the reverse proposal might equally 
be true but its proof is less obvious 0. 

The steps of the demonstration are as follows: 

• Choice of a dimer basis. Dimer decompositions are a bit uneasy because dimer 
coverings are not orthogonal and the entire family of dimer coverings is over- 



4 Miscellaneous remarks: i) a dimerized state does not fulfill property (A) because it has dimer-dimer 
long-range order (4-point correlation function), ii) the Spin-Liquid observed on the kagome lattice might 
not obey assumption A but it has most probably gapless singlet excitations. 
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complete. So we must suppose that in a first time we have extracted a non 
orthogonal basis of independent dimer coverings called generically \C). This 
basis should respect the translational invariance of the problem: which means 
that if |C) is a basis vector, so is T \ C) (where T is any unit step translation) 

Decomposition of the translation invariant ground-state. We have understood 
previously that our problem is a boundary problem. Let us draw a cut A 
encircling the torus created by periodic boundary conditions (see Fig. ||). This 
hyper-surface of dimension d — 1 cuts bonds of the lattice but there is no sites 
sitting on it. The position of the cut is arbitrary; but we may decide in order 
to follow closely our previous discussion to put it between spin L — 1 and spin 
of each row of the lattice. In the decomposition of |?/>o) on the dimer basis, let 
us sort the coverings in two sets A + and A_ according to the parity DA of the 
number of bounds going across the cut. This leads to a formal decomposition 
of \4>o) in two parts: 

\M = |V>o + ) + |V>o~) (12) 



where the vectors 



ip ) belong respectively to the sets A±. 



Let us now consider a finite system with an odd number of rows along the direc- 
tion of the cut, and a one-step translation 7a that crosses the cut. If \C) belongs 
to A + , Ta |C) belongs to A_ and reversely (Property (B)). The translational 
invariance of |^o) thus implies that both ipo^j and i/Jq^ are simultaneously non 
zero. With the assumption (A), it is easy to show that the two components 
are orthogonal. 

Excited state. One can thus build the state 

|^l,A> = |^o + ) " (13) 

Using the fact that | V^o) h as momentum zero, Property (B) implies that \ipi a) 
has a momentum kj\ in the direction of 7a- It seems reasonable to think that 
this property demonstrated to be true for samples with an odd number of rows 
is valid in the thermodynamic limit for any samples (this is clearly indicated 
by small size exact diagonalizations, see Fig. || for the 6x6 sample). Remark: 
I^IjA) is obtained from |^o) by changing the sign of the dimers crossing the 
boundary: in the case of a short range RVB state it is exactly the kind of local 
perturbation we were searching for in the previous paragraph. 

The demonstration is achieved by proving that IV'i.a) an d IV'o) have the same 
short range correlations and thus the same energy in the thermodynamic limit. 
More precisely only the phases of the long loops f\ in the expression of the 
energy (ipolTi-lipo) change sign and thus the variational energy of |V>i,a) only 
differs from the energy of |V>o) by terms of the order of 0(exp(—L/^)) where £ is 
the characteristic length of the loops. Such a construction can be done for any 
main direction of the lattice, which proves that the degeneracy of the ground- 
state in the thermodynamic limit is 2 d (a result already obtained by Read and 



'Loops with non-zero winding number around the torus. 
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Chakraborty, without reference to the wave vectors of the quasi degenerate 
ground-states) 

This completes our assertion that the LSM A conjecture is indeed valid in a very 
large number of situations whatever the dimension of the lattice. 



5 Miscellaneous remarks 

• We thus arrive at the conclusion that in the absence of long range order there 
is, strictly speaking, symmetry breaking of one step translations. Long range 
order in any of these antiferromagnetic systems implies symmetry breaking, 
BUT the reverse is false, the symmetry breaking described here does not imply 
long range order in a local order parameter . This property is of topological 
origin: the observation of this symmetry breaking would need sensitivity to a 
global observable or to boundary conditions (edge states). 

• It should be noticed that the results on the 6x6 sample seem to indicate 
that the conserved symmetry of the ground-state is in fact TS, where £ is the 
reflexion through a plane containing k^. 
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Figure 1: First eigenstates of the multiple-spin exchange model on a 6 x 6 sites sam- 
ple (ref. [fL0|1 , parameters J 2 = — 2 and J 4 = 1; the system is in a Spin-Liquid phase). 
Eigenstates with total spin S — and S — 1 are displayed. The symbols represent the 
spatial quantum numbers ( IZg, and cr = 1 are the phase factors taken by the many-body 
wave function in a rotation of the lattice of an angle 9 or in a reflexion symmetry). The 
ground-state (E = —142.867), belongs to the trivial representation of the space group: 
k = 0, 7Z2n = 1, H n = 1, cr = 1. The first excited states have wave-vectors (three-fold 
degeneracy). The finite size scaling strongly indicates that these states collapse to the 
absolute ground-state in the thermodynamic limit. The third and fourth eigen-levels in 
the S = spin sector do probably not collapse to the absolute ground-state fp7o| . They 



may be degenerate in the thermodynamic limit (4-fold degeneracy) and describe an S = 
bound-state just below the continuum of triplet excitations. 
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MSE Triangular J 2 =-2 J 4 = 1 




0.5 1 1.5 2 

0/(2tt) 



Figure 2: Low-lying spectrum of the 4x5 sample of the multiple-spin exchange model 
(same parameters as in Fig. [I]), as a function of the twist <fi. Squares are k = states in 
the Bfj, frame and circles stand for states with momentum k^ 1 (momentum n in the even 
direction). 




0.005 0.01 0.015 



(0/(2tt))4 

Figure 3: Ground-state energy as a function of 4 in the vicinity of <fi = (same parameters 
as in Fig. Q). The energy is not quadratic but rather proportional to 4 (vanishingly small 
stiffness). The line is a guide to the eye. 
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Figure 4: 2-torus with one cut A. 
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